Preprint typeset in JHEP style. - PAPER VERSION 



hep-ph/0010115 



Branching ratios for B^ s — > Jj^r} and B^ s —> r]£ + £~, 
extracting 7 from B^ s — > J/ip r\ , and possibilities for 
constraining Ciqa in semileptonic B decays 



P. Z. Skands 

HEP department - Niels Bohr Institute - Blegdamsvej 17 - DK-2100 Copenhagen 
zeiler@nbi.dk 



Abstract: Estimates of the branching ratios for B^g — > J/iprj and Bd s — » r/£ + £~ are 
obtained by SU(3) relation to B& — > J/ip K and B& — ► K£ + £~, respectively, as func- 
tions of the T)\ — r]s mixing parameter, Op. Based on these estimates, a discussion of the 
prospects for HERA-B, CDF-II, and ATLAS on these processes is given. The CP violation 
in Bd.s —> J/ipTj is analyzed in depth and a method to extract the angle 7 of the unitarity 
triangle is discussed. Finally, a possible method to constrain the Wilson coefficient C±qa 
from measurements on semileptonic B decays such as Bj — > K£ + £~ and B s — > 7/ £ + £~ is 
proposed along with a discussion of the prospects for future experiments and form factor 
calculations to reach the precision required for this method to be interesting 
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1. Introduction 



One of the most important processes relevant to the study of CP violation (see [1] for an 
excellent review of this topic) is the well-known "gold-plated mode" Bd — ► J/ipKs- Besides 
its usefulness in extracting sin(2/3) (see e.g. [2]), it has relatively recently been argued that 
this process can also be used to extract the angle 7 when combined with its ?7-spin partner 
B S ^J/^K S [3]. 

In this paper, an alternative to the gold-plated mode for the extraction of 7 is presented: 
Bd,s J/ip V • Due to the inherent problems in adapting factorization to non-leptonic 
decays we use the completely general and model independent method of quark topologies 
to analyze the structure of B^ s ~~ ^ J V > obtaining a relation between the Bd and B s 
amplitudes through SU(3) flavour symmetry (the approximate symmetry of u,d,s). The 
structure of the two processes is such that the phase, e* 7 , is CKM suppressed in B s — > J/iprf 
relative to Bd — > J/tpr) , and so effects of CP violation will be more easily visible in the 
Bd — ► J/ip V decay. However, the extraction of 7 from this process is plagued by the 
appearance of a normalization factor which cannot be determined directly. Through U 
spin relation, the CP averaged rates of B s — > J/ipV an d Bd — > J/tpr] can be combined to 
fix this normalization. Thus, 7 can be determined with a theoretical uncertainty depending 
only on SU(3) breaking corrections and the r\\ — rjs mixing angle, Op. Furthermore, the 
Bd amplitude itself is suppressed relative to the B s amplitude, the inverse of the situation 
in the Bd )S J/fp K$ case. The difference in production rates of Bd and B s mesons in 
the experiment compensates to some extent for this difference in the case Bd tS J/ipf) 
whereas it worsens the situation for Bd <s — > J/^Ks- 

As an addendum, a method to constrain C\qa is presented. In semileptonic decays, 
Cioa describes the effective coupling of the axial OPE operator Oiqa [1]. In models of 
new physics, its value generally deviates from the SM prediction, due to virtual particle 
contributions from New Physics particles present in such models. In this paper, a method 
is proposed which allows the elimination of large hadronic uncertainties caused by the 
presence of intermediate ip resonances by measuring distributions in semileptonic decays 
and using relations between them. It is, however, doubtful whether this method can 
find immediate application due to the high precision required both for theoretical and 
experimental input. 

The outline of the paper is as follows: In section 2, we use the method of quark 
topologies to analyze the CKM structure of the contributions to the Bd tS J/ipf) decay 
amplitude (see [4] for a recent update on this method). Noting the Zweig and SU(3) 
suppressions of certain topologies and assuming SU(3) symmetry of the strong interaction 
dynamics, a simple relation can be obtained between the Bd — > J ftp 7] and B s — > J/ip 77 
amplitudes. In section 3, the procedure proposed in [3] for extracting 7 from Bd :S — ► J/ipKg 
is adapted for Bd <s J/ip V , and the amplitude relation obtained in the previous section 
is used to obtain the normalization of the CP averaged Bd — > J/tpr/ rate. 

In Section 4, a simplified picture of the Bd tS J/ip V process is adopted and SU(3) 
symmetry is invoked to obtain a relation between the amplitudes of Bd — > J/tprj , B s — > 
J/iprf, B s — > J/ip Ks to the already measured Bd — > J/ipKs- In section 5, essentially the 
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same is employed to obtain Bd, s — ► V form factors from those for Bj — > K calculated by 
Light Cone Sum Rules (LCSR) in [5JI6] 1 . Using these, an estimate for the B^ s — > r]£ + £~ 
branching ratios can be obtained. Both the B^ g — > J/i/jt) and the B^ s ??^ + ^~branchings 
depend on the degree of octet-singlet mixing in the r\ system, expressed through the mixing 
angle, Op. 

In section 6, a simple method is proposed to constrain C\oa, the axial semileptonic 
Wilson coefficient. We replace the theoretically poorly known quantity Cgy in the B& — ► 
Kt + t~ amplitude by measurable decay distributions for the B^ — ► Kfi + fi~ process, yield- 
ing Cioa as a function of the total branching to Ktt, the differential branching to K/j,/j,, 
\V t * s V t b\ 2 , and the form factors /+ and /_. Estimates in various SUGRA models and the 2 
Higgs doublet model are considered, and the application of the same procedure to the case 
B s — ► i] is considered. Concluding remarks and outlook are given in section 7. 



2. Analysis of B d:S — > J/ipr) 



-ftO 



The time-independent transition amplitudes for B° and B" states into the final CP eigen- 
state, fcp, are parametrized as in [3]: 



{f CP \H\B°)=A 



1 - a q e t0 "e^ 



= N q z q 



(f CP \H\B q )^A q = V N q 



l-a q e w "e-^ 



= T]N q Z q 



(2.1) 
(2.2) 



where r\ is the CP eigenvalue of fcp, is a strong phase, and q G {d, s}. The amplitude 
at time t for an initial B/B meson to decay then becomes: 



\A q {t)f 



\A q (t)? 



\N n 



[{R L + e pq B} L )e- TLt + (R H + e pq R l H )e- THt 



+2e 
liVJ 2 



-Tt 



{A D + e pq A D ) cos(AMt) + (A M (1 - -f-) sin(AMt) 



(2.3) 



(R L + e pq R l L )e~ VLt + (R H + e pq R l H )e- THt 



-si 



-2e 



-Tt 



{A D + e pq A D ) cos(AMt) + (A M {1 + ^) sin(AMt) 



(2.4) 



where we define rate functions as in [3] , except that we here give the formulae to first order 
in the small parameter e pq = ||| 2 — 1, where p and q are the standard factors parametrizing 
the B meson mass eigenstates in terms of flavour states [1]. As there is, however, only small 
possibility for new physics to result in a large e pq we henceforth ignore this parameter. The 



1 For a more recent calculation, see [7]. 
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rate functions entering the above are denned by (4> q is the B q — B q mixing phase): 



Rl 


= 


-(\z q \ 2 + \z q \ 2 + 2r ] Re{z* q z q e^}) 


(2.5) 


Rh 


= 


l(\ Zq \ 2 + \z q \ 2 -2 V Re{z* g z q e^}) 


(2.6) 


A d 


= 




(2.7) 


A 

A M 




—7]lm{z q z q e Vq ) 


(2.8) 


Rl 

Li 




--(\z a \ 2 + ■nReUnZae i4 ' q }) = --(R L + An) 


(2.9) 


R\j 




--(\z q \ 2 - rjRB{z*z q e^}) = --(R H - A D ) 


(2.10) 


M> 




2 2 \2 / 


(2.11) 






i(|z g | 2 + r/Re-^e^}) = ^ + ^) 


(2.12) 


R H 




|(l^| 2 " rjRB{zp q e^}) = \(R H - A D ) 


(2.13) 






\\^? = \(\{Rh + Rl)A d \ 


(2.14) 



With regard to the final state itself, a few comments are necessary regarding the octet- 
singlet mixing in the rj system. This mixing, parametrized by the mixing angle Op, is still a 
controversial issue (for current experimental values, see e.g. [8]), so rather than using some 
specific value for Op, we rewrite the r\ wavefunction in the following way: 

rj = —={uu + dd — 2ss) cos(0p) -j={uu + dd + ss) sin(0p) 

V6 V3 

= N v (uu + dd) + S v (ss) (2.15) 
_ cos(flp) sin(gp) _ -sin(flp) 2cos(g P ) 

N " = ~ve 7T ^ = —3 W (2 - 16) 

which is the definition we shall use in the following. 



2.1 The time-independent amplitudes 

The expressions given above for |^4 g (t)| 2 and |^4 g (t)| 2 depend on the time-independent am- 
plitudes parametrized by eqs. (2.1) and (2.2). In non-leptonic processes, these amplitudes 
can generally not be evaluated by any method relying on the factorization approach due 
to final state interaction (FSI) effects. In the present work, however, we do not need an 
explicit calculation. Rather, we wish to obtain a parametrization of the B& and B s am- 
plitudes that will allow us to relate them by SU(3) flavour symmetry. To arrive at such 
a parametrization, it is sufficient to use the method of quark topologies which does not 
require the ability to solve the full theory and which allows a systematic classification 
of long-distance contributions (for a recent update on this method, see [4]). With the 
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PENGUIN EMISSION DOUBLE PENGUIN ANNIHILATION 

Figure 1: Quark topologies in _B djS — > J/ip-q. The gray blobs denote the contracted short-distance 
parts. 



topologies listed in figure 1, we obtain the following 2 : 

A(B d - J^r? ) = N v {\ c bd A d + A&5 d ) + \ u bd N v 
A{B S -^J/ 1 p v ) = S v (X c bs A s + XtB s ) + A^iV, 

with the definitions: 



A bd 



D d + T^rCd + £d 
\d 

Ds + ^Cs + Cs 
A bs 



A x 



B x = 
D x = 

Cx = 



<QST> re + <e^> E + E<^>pI 

\«1,2 / E A2 

((Q^f) EA1 + (Q^f) DPA + «2&*> E 

1 \^pen I fa? \^pen I fa? iVn \~ c pen I 



^(Qpera ) DPA ) 
or / 



-pen / E A2 



1 \^pen / F a? \~°pen i fa7 iVr? \ -°pen I 



pen / EA2 



(2.17) 
(2.18) 

(2.19) 

(2.20) 
(2.21) 

(2.22) 

(2.23) 



where A^ = V*, b V q / q , (Qf ) x denotes the insertion of the OPE operator Qi having external 
quark lines buud into topology x, and (Q. 



(u)d x 
pen 



denotes the combined contribution of the 



QCD penguin operators 63-6- The quark content is denoted (q)d,(q)s where q is the 
flavour of the qq pair coming from the gluon and d or s are from the flavour-changing 
6-quark transition. It should be mentioned that, relative to the current-current operators 
Qi,2, the electroweak (QCD) penguins contain an extra power of olem («s)- In neither 



The computational details can be found in an unpublished project. Please contact the author if a copy 



is needed. 
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the Bd nor the B s case are the dominant current-current contributions CKM suppressed 
relative to the penguins, and so we expect |4ewp|/|A7c| = 0(1O -2 ) = |Aewp|/|^4qgtdp|- 
I have therefore neglected elelectroweak penguin contributions in the above analysis. 

Noting that the £ and £ terms are SU(3) suppressed and that the D terms are OZI 
suppressed (see figure 1), we neglect these and obtain: 



A d = A(B d 
A s = A(B S 



j/V ?? ) = Nr, (v; b v cd A d + v: b v ud B d ) 

= N d [l - aae^e^ 

j/i>v) = s v (v; b v cs A s + V* h V us B s ) 



= N„ 



A 2 JOs Jl 



(2.24) 
(2.25) 



with 



N d = -N r ,AX 3 A d , a d = 
N s = S ri AX 2 A s fl-y | . a H = 



2 



R b (l- 
R b (l - f) 



and the Wolfenstein parameters [3]: 

A ee \V US \ = 0.22, A ee ^|y c6 | = 0.81 ± 0.06, 



Bd 

Ms 

As 



Rh = 



d = Aig{§*} (2.26) 
9 S ^-Arg{^} (2.27) 



A 



V ub 



0.41 ±0.07(2.28) 



With eqs. (2.24) and (2.25) we have recovered exactly the form of eq. (2.1) by which we 
parametrized the time-independent amplitudes, but in a form that explicitly separates the 
CKM structure from the strong amplitudes. This is of essential use below. 

3. Extracting 7 from B djS — > J/ipr] 

We here adapt the method proposed in [3] for the decay considered here. Defining the 
time-dependent CP asymmetry by: 



acp 



\A q {t)\* - \A q (t)? 



\A q {t)\ 2 + \A q (t)\ 2 
and inserting the above expressions, one obtains [3]: 

Afp cos (AM t) + A^f sin(AM t) 



acp{t) = 2e 



-vt 



e -T H t + e -T L t + j{ AT ( e -T H t _ e -T L t^ 



with 



\dir 



a mix 



2A 



D 



2d q sin 7 sin 9 q 



CP = 



Aat = 



Rh + Rl 

2A m 
Rh + Rl 

Rh-Rl 
Rh + Rl 



1 — 2d q cos 7 cos 9 q + a q 
sin (f) — 2d q sin(7 + (f> 

77 



sin (p — 2d q sin(7 + (f) cos 9 q + a q sin(27 + 4>) 

1 — 2d q COS 7 COS 9 q + 0% 

cos 4> q — 2d q cos(7 + (f) q ) cos 9 q + d q cos(27 + 



-v- 



1 — 2d q cos 7 cos 9 q + d\ 



n2 



(3.1) 

(3.2) 

(3.3) 
(3.4) 
(3.5) 
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a d 

e s + i80 c 



; for B d 
: for B s 



J/lfjT] 

J/ifjr) 



; for B d ->■ J/tpi] 
; for B s — ► J/i^rj 



(3.6) 
(3.7) 



As a s is suppressed by y^p-, we expect a very small direct CP asymmetry in the B s 
process, thus we shall use the asymmetries in B d — > J/ip r] combined with R d and R s 
for extracting 7 (R q = \{R q H + R q L )). It should be mentioned that, due to the smallness 
of AI^, the 'observable' »4Ar will only be measurable for the B s system. This has no 
direct consequence on the analysis presented here. Only two of the three asymmetries are 
independent (see below), and so a measurement of A^p and A™p is sufficient. 



{A d S r P f + {Acpf + {A^f = 1 



(3.8) 



which has been checked also to be valid when going to e pq / 0. 

The observables A^fp and A™p do not depend on the normalization, |A^| 2 , and can 
in principle be obtained by fitting to the CP asymmetry. This yields two equations (3.3)- 
(3.4) in the three "unknowns": a d , 6 d , and 7 (taking the mixing angle, 2(3, to be known 
beforehand). Thus, we need one more observable. Measuring the CP averaged rate yields 



<r,> = n 2 x |i\g 2 x R q 



(3.9) 



where II2 is the 2-body phase space, and the normalization factors for the Bd and B s modes 
are given by eqs. (2.26) and (2.27). A priori, we cannot determine this normalization, but 
assuming now that the strong interaction dynamics is the same for the B d and B s modes 
([/-spin symmetry), we have a d = a s = a, 9 d = 6 S = 9, and taking kinematics into account: 



Ns 
N d 



X 2 



A 2 



A d 



A 2 



A 2 A(m| d ,m 2 / ^ ) ,m 2 ) 



(3.10) 



with the A function being the standard one used in kinematics. Thus, forming the ratio 
R s /R d = (R s H +R s L )/(Rj I +Rf i ), the strong dynamics cancels out, and we find by combining 
eq. (3.9) with R = \{Rh + Rl) = \ (1 — 2a q cos 7 cos 1 



&d 



1 - 2j^a s cos7cos# s + 

1 — 2a d cos 7 cos 6 d + a 2 d 

(T s ) \N d \ 2 m| s / A(m 2 



2^2 



, m 



, mz 



(T d ) \N S \ 2 m% d \\{m 2 BsJ m 2 J/r m 2 
x2 (T s )\N v \ 2 m 3 Bs (\(m 2 



B d ' m J/i, 



1/2 



to 2 )' 



(T d ) |^| 2 m| d y\{m 2 Bs ,m 2 J/r m 2 ) 



3/2 



(3.11) 



This relation furnishes us with the last observable we need for determining a, 9, and 7 as 
functions of the mixing phase, <fi. In the Standard Model, this angle is negligibly small for 
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6 . S > B 

s d 




d 



Figure 2: A simplified picture of the four SU(3) related decays Bd, s — » <//V' r 7 an d — * J/ipK°. 

the 5 S — £> s system whereas it is 2/3 for the .B^ — B d system. Thus, the extraction of 7 
requires (3 as an input parameter. 

To summarize: The two observables in eqs. (3.3) and (3.4) are to be determined for 
the decay B d — > J ftp rj by measuring the CP asymmetry and fitting to the time-dependent 
decay amplitudes. This will require tagging and will fix a contour in the 7 — a plane. The 
last observable is provided by H = R s /R d which does not require tagging as it depends 
only on the CP averaged rates. Together with e.g. A™™ , we fix another contour in the 
7 — a plane. The intersections of these two contours will fix both a and 7, to a theoretical 
precision depending on SU(3)-breaking corrections and the accuracies on the rj and B d 
mixing angles, 6p and 2/3. 

As we have chosen the same terminology and parametrizations as [3] the contour 
equations in that paper are directly applicable, and we do not repeat them here. 

4. Branchings for B diS — > J/tprj 

We now present a simplified picture of the B d jS — > J/V"7 transitions in order to obtain an 
estimate for the branching ratios. We again take the strong interaction dynamics symmetric 
under SU(3) transformations and further make the simplifying ansatz of fig. 2 that the 
insertions into emission topologies represent the dominant contributions, resulting in the 
parameters a d and a s of the previous sections being negligibly small (we really only have 
to make this crude assumption for the B d — ► J ftp rj and B s — ► J/tp Kg processes. In 
B s — > J '/ip 77 and B d — > J/ip K$ (see [3]), we can justify neglecting these terms because of 
the A 2 suppression). With these assumptions, the amplitudes in fig. 2 can differ only by 
CKM factors, kinematics and factors coming from the hadronic wavefunctions. Indeed, 
eqs. (2.24) and (2.25) become (separating kinematics and dynamics): 

A(B d - J/V 77 ) = N v X c bd A d = N v X c bd ( PBd + Pv re^F(B d - J/^r,) (4.1) 
A(B S - J/^r?) = S v X c bs A s = S v X c bs (PB a +p v Ye^F(B s - J/^r,) (4.2) 



7 



where p and e are momenta and polarization vectors, respectively, and the form factors, 
F, parametrize the strong dynamics. 

For the B ds — > J/ipKs I use the amplitudes in [3] with the same assumptions as above 
and with a slight modification of the author's notation. 

A(B d - J/^K S ) = \ c bs A d = X c bs (p Bd +p Ks Y^F(B d - JH>K S ) (4.3) 
A(B S - J/^ tf s ) = = A^Cpb. + PKsT^HBs - J/^ ^s) (4.4) 

With the assumption that the strong dynamics is SU(3) symmetric, the form factors cancel 
out when forming ratios, and we arrive at: 

A(B d ^J/i>7]) _ ^\ c bd (p Bd 



A(5 d J/^ K 5 ) " XI ( PBd + p**)^ 



A(B d - J/^ " (p Bd + s )"e„ 



(4.7) 



A(£ d - J/^ K 5 ) Ag s (p Bd + PX S )% 

where the \/2 comes from the translation from K° to Kg- In addition to the modes we are 
interested in, we have written down A(B S — > J/i/jKs) as a bonus. 

By using these relations, it is straightforward to obtain the branching ratios for B d s — > 
J/^ 77 and S s — ► ^Cs in terms of that for B d — ► J/^ Kg with the measured value 
Bi2(.B d -► J/ipK°) = 2BR{B d -► J/4>K S ) = (8.9 ± 1.2) x [8]. Inserting this value 

yields: 



2 2\ \ 3/2 

^,->/,„ -9X10 ^ (S^St) 



/, ,o m H / A(m?, ,m 2 7 , , .mi) \ ^ 2 

BRb^j,^ = 9 x 10- 4 |^| 2 ^f A( L£>^4) ) (49) 



BR Ba _+j/^ K o = 9 x 10" 



4 




2 H d / 









2 2 \ \ 3/2 

' l J/ip< m K) 



M m B d > m J/4> m K 



(4.10) 



As mentioned in section 2, there is some controversy as to the precise value of the rj mixing 
angle, Op, which determines N v and S v . Rather than adopting some specific value, we 
have varied the parameter between —20° < Op < —10°, producing the results shown in 
table 1. The #p-dependence between these limits is linear to a good approximation. The 
uncertainty on these branching ratios is roughly 40%, slightly more for B d — ► J/V"7 and 
B s ^ J/^K S . 

According to the PYTHIA simulation for HERA-B [9], the production rate of B d 
mesons is about five times greater than that of B s mesons, and so we will expect to see 
about 5 times more B s decays than B d decays in the experiment for Op = —20°. For 
Op = —10°, we expect about 10 times more B s decays than B d decays. This situation is 
slightly better than for the B d s — > J/ipKg decays [3] where we expect to see around 250 
B d events for each B s event, and so the statistical error on the B s events is going to have 
a larger influence on the precision with which 7 can be extracted. In both strategies, the 
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B d ^J/^ V 


B s ^j/ip v 


HERA-B 9 P = -10° 
(untagged, /yr) 6> P = -20° 


100 
125 


850 
560 


CDF II 6>p = -10° 
(untagged, 2fb _1 ) 6 P = -20° 


700 
1000 


5.8 x 10 3 
3.8 x 10 3 


ATLAS Bp = -10° 
(tagged, 30fb _1 ) P = -20° 


4 x 10 4 

5 x 10 4 


3.5 x 10 5 
2.5 x 10 5 


77 not reconstructed (— » factor 10% — 20%) 



Table 2: Estimated number of reconstructed Bd. s — » J/ip f] events at HERA-B, CDFII, and 
ATLAS. 



asymmetries are to be determined for the mode with least statistics, so a few more factors 
are definitely of use. Given the branching ratios, it has been estimated how many events 
will be seen by HERA-B, CDF-II, and ATLAS, based on their simulations for B d — > J/t^Ks 
[10] [11] [12]. Results are listed in table 2. As the r\ reconstruction efficiency is not known 
at present, the numbers presented here are without rj reconstruction included. A loose 
estimate of this efficiency is 10-20%. 





Bp = -10° 


Bp = -20° 


BR(B d ^J/i> V ) 
BR(B s ^J/^rj) 


1.1 x 10~ 5 
5.0 x 10~ 4 


1.5 x 10~ 5 
3.3 x 10~ 4 


BR{B S -► J/^ U ) 


5.28 x 10~ 5 



Taking the rj reconstruction efficiency 
into account, it is certain that HERA-B 
will not be able to access this mode (the 
number given is for the machine run- 
ning at full luminosity), it is an open 



Table 1: Estimated branching ratios for Bd. s — * question whether CDF-II will be able 
J/iprj and Bd — > J/tpK°. to get it (depending on how much lu- 

minosity they get before LHC, and whether 
they improve their trigger efficiency [11] [13], and it is certain that ATLAS will access it 
within the first three years of operation. 



5. Branchings for B^ s — > r]£ + £ 

The branchings for B s — > rj can be obtained by using that the process is related 

by the approximate SU(3) flavour symmetry to Bd — ► Kt r t~ whose spectrum has been 
calculated in [6]. Due to the close similarities between r\ and K mesons (pseudoscalars, 
similar masses), and going to the SU(3) symmetric limit, we merely need to replace B — > K 
by B — > rj form factors and to take CKM factors into account. At present, no reliable form 
factor calculations for B^^ ~^ V exist. 

In the following, we estimate the form factors for B s — > rj by SU(3) relation to the 
B — > K form factors presented in [6], effectively resulting in a multiplication of the form 
factors by S v (see eq. (2.15)). The calculation of Bd — ► rj form factors is essentially identical, 
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and so we omit an explicit calculation. 



(r)(p v ) \s^b\B s (p B )) = 

Sr,\f+(s)(p B +P V )» + f-(s)(pB- P r,)^} 

(ri(p v )\sa llu q u b\B s (p B )} = 

Sr,i[s(p B + P V )» ~ {m 2 B - m^q^ J^l^ 



(5.1) 



With these form factors, the total branching ratios for Bd tS — ¥ r\£ £~ as a function of Op has 
been calculated using eq. (6.1) integrated over the kinematical region without intermediate 
tp resonant states included in Cg^. 

A conservative theoretical uncertainty on 
this calculation is given by « 40% from 
SU(3) breaking corrections, and ~ 25% from 
uncertainties on the initial B — > K form 
factors, yielding a total uncertainty around 
50%. Results are listed in table 3 for the 
case i = \i. Based on simulations from the two experiments of B —>■ K*£ + £~, the num- 
ber of reconstructed B s — > rj fi + ' \l~ and Bd — > rj events have been estimated, again 
without taking the rj reconstruction efficiency into account (table 4). 





Op = -10° 


P = -20° 


B s -> 7]fl + fl 


3.0 x 10~ 7 


2.0 x 10~ 7 


Bd^V 


7.7 x 10~ 9 


1.0 x 10~ 8 



Table 3: Branching ratios for Bd. s — > r/£ + £ 



CDF-II 

(2 fb^ 1 ) 



9 P 



-10° 
-20° 



ATLAS e P = -10° 
(30 ft" 1 ) e P = -20° 



B d 



10 



20 
25 



60 
40 



160 
105 



rj not reconstructed (— > factor 10% — 20%) 



Table 4: Estimated number of reconstructed Bd, s — * rj£ + £ events at CDF-II and ATLAS. 



6. Constraining C i0 A in semileptonic B decays 

The differential branching ratios of Bd s — ► K£ + £~ and B s d —> rj£ + £~ can be written as 
[6]: 



dT GWm 5 



ds 



B IT/* 



2 10 7T 5 



V t l td V tb \ z u(s) 



(\A'\ 2 + |C'| 2 )(A 



U(S 



+\C'\ 2 4mj(2 + 2m£ - a) + Re(C"l>' )8m|(l - m£-) 
+ |,D'| 2 4mp] 



(6.1) 



with 



X K}V = X(l,m 2 KjV ,s) u(s) = dX(l-4-^) 



(6.2) 



10 



and the coefficients A', C', D' given in terms of Wilson coefficients and form factors as [6]: 

A'(s) = Cf(s)f^(s) + - 2 ™ b Cff^(s) (6.3) 

m B + m K ^ 

C'(s) = CW+ •"(*), D'(s) = CioA/5"(a). (6.4) 

Variables with " have been normalized by m 2 B , A is the standard kinematical function, and 
s is the energy in the CM of the lepton pair. The Cj entering these expressions are the 
Wilson coefficients of the OPE operators contributing to semileptonic B decays: 

= ^m b s a a^{\ + 75 )b a F^ u Oiy = (s a b a ) v . A (U) v O e 10A = (s a b a ) v _ A (U) A 

(6.5) 

Cf^ = C-j + C5/3 — Cq and contains corrections due to intermediate quark loops in 
the decay (see e.g. [6]). 

For zero lepton mass [l = e, fj) only the ( | v4 r | 2 -h | C r | 2 ) term in (6.1) remains. Measuring 
dT /ds for e.g. the B — > K/i + /i~ thus yields a way of experimentally determining this sum. 
This is attractive since A' is excessively difficult to calculate theoretically, due to large 
hadronic uncertainties caused by the existence of intermediate ip resonances entering C^. 

In contrast to this, the terms in (6.1) containing C and D' depend only on C\q A and 
form factors. Eliminating A 1 from the r distribution by isolating the term (| A'\ 2 + |C"| 2 ) in 
the muon distribution and inserting it in eq. (6.1) for t = r yields the following relation 3 , 
where we have integrated over the kinematical region of the r spectrum: 



r{m B -m K ) dB r{m B -rh K ) I ^2 dB ( 2rhl\ IA 

/ ^ ds= A/ 1 " - ^ ^ (l + —^)ds 

J Am 2 as y 4A 2 V s ds V s ) 

Ct 2 n 2 m b f(m B -rh K ) 2 

MC W A\ 2 ^^\V*V tb \ 2 T B / F(s)ds (6.6) 

z " JAm 2 

where tb is the B meson lifetime and 



F(s) = - {f 2 4m 2 T (2 + 2m 2 K - s) +f + f^m 2 T (l - m 2 K ) + fHm 2 s) 

To illustrate the order of magnitude of the terms entering eq. (6.6), we give their SM values 
below (with the form factors presented in [6]). 

/ = 1.35 x i(T 7 (6.7) 

J Am 2 

r(m B -m K ) 2 jg 

/ ds ^l-Am 2 /s (l + 2m 2 /s) = 0.90 x 10~ 7 (6.8) 

J Am 2 . ^S 

2-5 r(m B -m K ) 2 

ICWI 2 010 5 B \V t * a V tb \ 2 TB / ds F(s) = 0.45 x lO" 7 (6.9) 

^ ft J Am 2 



3 The result for B s — ► rj is completely analogous. 
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Due to the close similarity between K and 77, these results are not changed substantially 
when considering B s — > 77 Using very precise form factor calculations, F(s) can 

be evaluated theoretically to within perhaps ±20% uncertainty. At second generation 
machines, it is reasonable to expect measurements of BR(B — ► Kt + t~) and the muon 
spectrum to about 10% accuracy. Isolating C\qa in the above formula would then yield 
an overall uncertainty of approximately ±20%. In the SUGRA models investigated in [14] 
(minimal and non-minimal, tan/3 = 2 and tan/3 = 30), C\qa lies within ±10% of the SM 
value, and so a distinction is not yet possible for these cases. In the 2HDM 4 , C\qa receives 
contributions from charged Higgs bosons, and we have used [15] for values of the Yukawa 
coupling of the charged Higgs to the top quark, < X tt < 0.3 and charged Higgs masses 0.2 
TeV < Mjj± < 1 TeV. The case of maximal deviation (Mjj± = 1 TeV, Xu = 0.3) produces 
only about 5% deviation. Thus we see that some of the most common models cannot 
be ruled out by a measurement such as the one proposed here in the near future. Both 
form factor calculations and experimental precision must be improved before this method 
becomes viable. Again, this result is not substantially altered by having a Bs initial state 
and an 77 in the final state. 

7. Conclusion 

Using the method of quark topologies and invoking SU(3)-suppression of certain topologies 
as well as SU(3) symmetry of the strong interaction dynamics, it is possible to obtain a 
parametrization of the Bd )S J/ipr) amplitude which allows the extraction of the angle 
7 of the Unitarity Triangle through measurement of CP asymmetries in the B^ — ► J/tprj 
process and measurement of the CP averaged widths of the B^ s — ► J/ip rj processes 
together. Electroweak penguins cannot lead to any problems in this extraction. The 
measurement of asymmetries require tagging whereas the CP averaged rates do not. 

Additionally, an estimate of the branching ratios for B& — > J/iprj has been presented, 
yielding an expected branching ratio of around 1 x 10 -5 for B^ —> J/iprj and around 5 x 10~ 4 
for B s — > J/ipr] (taking 9p ~ —20°). Due to the difference in production ratios of B^ and 
B s , this means a factor 5-10 more B s decays than B^ decays in the experiment, depending 
on the value of 9p. This is an improved situation relative to the B^ s — > J/ip K$ decays, 
where a factor of around 250 difference is expected. It has been estimated that HERA-B 
will probably not be able to see Bj — > J/ip rj decays whereas the Tevatron experiments 
CDF-II/D0 perhaps stand a chance of extracting of 7 from B^ s J/ipr) ■ Finally, the 
ATLAS experiment should have sufficient statistics to allow the extraction of 7 even taking 
the most pessimistic approach to the uncertainties in the estimates of branching ratios and 
production rates presented in this paper. 

As an addendum, a strategy to constrain C\qa in semileptonic B decays with small 
hadronic uncertainties has been proposed. It has been found to require both high theoretical 
and experimental precision, of a kind that will not be available in the near future. 

4 Two Higgs Doublet Model 
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